Three-Dimensional Light Bullets in Arrays of Waveguides 
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We report the first experimental observation of 3D-LBs, excited by femtosecond pulses in a system 
featuring quasi-instantaneous cubic nonlinearity and a periodic, transversally-modulated refractive 
index. Stringent evidence of the excitation of LBs is based on time-gated images and spectra which 
perfectly match our numerical simulations. Furthermore, we reveal a novel evolution mechanism 
forcing the LBs to follow varying dispersion/diffraction conditions, until they leave their existence 
range and decay. 
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Since their theoretical prediction (Tj, light bullets (LB) 
constitute a frontier in nonlinear science. These solitary 
waves are spatiotemporally localized and, in particular, 
find their manifestation in three-dimensional (3D) sys- 
tems governed by the Nonlinear Schrodinger Equation 
(NLSE), such as quantum gases [2|, evaporating black- 
holes [3|], and solid state physics In contrast to the 
one-dimensional case (e.g. optical fibers (Hf), where the 
NLSE is fully integrable and supports soliton solutions, 
higher dimensional solitary waves and especially LBs are 
not solitons in the strict sense of integrability of the dy- 
namic equations, thus being subject to instability [6| and 
limited exsitence range Their appeal as particle-like 
wavepackets triggered a two-decades-long research for a 
stabilization mechanism enabling full-dimensional, non- 
linear light localization. Theory shows that LBs can 
be stabilized by a variety of experimentally motivated 
modifications of the NLSE, such as saturation of the 
nonlinearity [8], higher order diffraction/dispersion [9], 
or nonlocal nonlinearity [10]. Transversally modulated, 
nonlinear media [11], e.g. arrays of evanescently cou- 
pled waveguides, have also been predicted to support 
stable LBs |3, EH- Despite theoretical progress, LBs 
eluded experimental observation. Experiments designed 
to observe spatiotemporal localization in nonlinear pla- 
nar media [13] revealed spatiotemporal compression, but 
the complexity of the observations could not be read- 
ily associated to LBs. The complexity is due to exist- 
ing optical media which support LBs only for conditions 
where effects beyond the Kerr-nonlinearity are influen- 
tial, e.g. self-steepening and intrapulse Raman scatter- 
ing. For this reasons, spatio-temporal solitary waves were 
first observed in x^^^ media where higher order effects 
were made negligible by artificially enhancing temporal 
dispersion by means of the tilted pulse technique Q. 
We remark that even in far-from-ideal systems the con- 
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FIG. 1. (a) Layout of the experiment. A focused 170 fs in- 
frared pulse centered at Ao = 1550 nm excites the central 
waveguide of a sample of the waveguide array (b). The out- 
put radiation of the sample is characterized by an infrared 
camera (IR-Camera) and by a spatially resolved optical gat- 
ing via frequency mixing in a 25-/im-thin BBO crystal with a 
delayed 60 fs probe pulse at Ap = 800 nm. (c) Experimental 
discrete linear diffraction pattern at the end of a 50 mm-long 
sample. Parameters of the sample: array made of 91 silica 
cores embedded in a Fluoride-doped silica glass; lattice pe- 
riod oi d — 33.2 ± 0.4/im, a core radius of r = 9.8 it 0.2/im, 
and an index contrast of An = 1.2 • 10~^. 



cept of solitary waves is a useful tool, allowing the un- 
derstanding of complex nonlinear phenomena such as op- 
tical rogue waves [l5|, and super-continuum generation 
m, [13. In fact, the name 'quasi-soliton' was recently 
used to describe nearly stationary spatiotemporal wave- 
packets propagating in planar waveguides arrays [18]. 

In this Letter, we report the first observation of 3D-LBs 
in a two-dimensional array of coupled waveguides. We 
have found that due to higher order effects, LBs evolve 
following varying dispersion/diffraction conditions until 
they leave their existence range [7] and decay. 
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FIG. 2. Experimental observation and simulation of LBs for a 40-mm-long waveguide array. Spatial infrared camera image 
(first column). Temporal cross-correlation trace of the central waveguide (second column) with the original experimental data 
(blue line) and after deconvolution of the probe pulse profile (red line) . Experimental spatiotemporal images (third column) and 
corresponding simulated spatiotemporal images derived (fourth column). Surfaces in the spatiotemporal images are normalized 
isointensity levels as indicated in the legend in (d).The three rows correspond to the three different input peak powers. 



The experimentally investigated system consists of an 
hexagonal array of evanescently coupled single mode fi- 
bres, where a focused femtosecond pulse is used to excite 
the central waveguide. A microscopic image of the array 
cross-section is shown in Fig. 1 along with a sketch of 
the experimental setup. The high regularity and symme- 
try of the observed diffraction pattern (Fig. 1(c)) proves 
the unprecedented quality of the array [l9|, and is a key 
for the observation of LBs, which was possible only with 
a high-resolution, spatially-resolved optical pulse cross- 
correlation techniques [20, 21] (see Fig. 1(a)). 

Temporal dispersion is anomalous around the excita- 
tion wavelength Aq = 1550 nm, but terms beyond the 
parabolic approximation cannot be neglected. Spatial 
discrete diffraction originating from inter- waveguide cou- 
pling [ll| is also inherently wavelength dependent. 

Figure 2 summarizes the results of experiments where 
the central waveguide of a 40 mm-long array of fibres 
was excited with 170 fs pulses. The spatial profile, the 
intensity cross-correlation trace of the central waveguide, 
the experimental, and simulated spatiotemporal profiles 
of the output wavepackets are shown for three different 
input power levels. At low input power (Pq = 0.2 MW) 
light spreads into the neighbouring waveguides while the 
pulse profile broadens (Fig. [2]a-d). By raising the input 
power to Pq = 0.4 MW, we observe a sharp localization of 



light in the initially excited waveguide with a pulse com- 
pression below the resolution of the cross-correlator (^ 60 
fs, see Fig. [2]e-h), while visible radiation is observed even 
by the naked eye. For higher powers (Pq = 1.0 MW), the 
light localization is enhanced and multiple short peaks 
are observed in the cross-correlations, suggesting the for- 
mation of bullet-like entities. 

Numerical simulations are based on the Unidirectional 
Maxwell Equations [16], which describe accurately the 
propagation of light in the array. Dispersion to all orders, 
wavelength-dependent discrete diffraction, Kerr nonlin- 
earity, intrapulse Raman scattering and self-steepening 
are accounted for. The remarkable agreement with the 
experiment (compare Fig. 2 (d), (h), and (n) with Fig. 
2 (c), (g), and (m)) did not require the adjustment of 
free parameters, thus justifying the use of numerical sim- 
ulations to show that the observed localized wavepackets 
are indeed LBs. Fig. 3(a) displays the evolution of the 
temporal profile of the light in the central waveguide of 
the array for an input power of Pq = 0.9 MW. Within the 
first 15 mm of propagation, the excited pulse gets spa- 
tiotemporally focused and splits into several fragments 
with (FWHM) durations ranging from 15 to 30 fs. The 
central wavelength of the fragments increases rapidly to 
larger values in proportion to their initial intensity, thus 
accounting for the observed distribution of group veloc- 
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FIG. 3. Numerical propagation of LBs in the central waveg- 
uide for Po = 0.9 MW. (a) Wavepacket evolution in the time- 
propagation plane. The diffraction and dispersion length cor- 
responding to a tfwhm = 25 fs pulse for the minimal and 
maximal central wavelength of the LB are overlaid to the 
graph for comparison, (b) Plot of pulse duration and central 
wavelength for the brightest LB seen in (a). 



ities. The brightest pulse experiences several collisions 
and eventually propagates as virtually stationary LBs to 
a distance of 2: = 40 mm (see Fig. 3(b)), beyond which 
the wavepacket suddenly spreads in space and time. Be- 
cause the central wavelength of the LB is not constant 
and the characteristic dispersion and diffraction lengths 
(^disp, ^diff) are wavelength dependent, we can gauge 
the roboustness of the localization of light in the array 
by defining an average dispersion and diffraction length 
(^disp) and (Ldiff) for the brightest pulse (see definition 
at the end of the Letter). According to this definition, 
nearly stationary propagation is achieved over 1.9 (i^diff) 
and 9.0 (i^disp)- We also proved that the light localization 
is not due to conical waves [22] by numerically verifying 
that the wavepackets do not self-heal. 

This scenario of the LB evolution has been confirmed 
experimentally by means of Frequency Resolved Optical 
Gating (FROG) [2lj obtained from a set of spatiotem- 
poral images, each recorded by inserting after the sum- 
frequency crystal a 10-nm interferential filter chosen form 
a set of 6, with central transmission wavelength ranging 
from 510 to 560 nm. This is equivalent to sampling the 
time-frequency plane of the LBs between A = 1423 nm 
and A = 1870 nm. In the left column of Fig. 4, we 
show the experimental FROG traces for the pulse in the 
central waveguide of three samples of the array (length 
L = 25, 40 and 60 mm, Pq = 0.9 MW). Frequency- wise 
interpolation was used to display the data. Two LBs 
with central wavelengths of 1650 nm and 1810 nm are 
observed at 2: = 25 mm. Only a LB with Ac = 1790 nm 
is observed at 2: = 40 mm. The experimental data can be 
interpreted by the corresponding simulated FROG traces 
(Fig. 4, right column), which accurately reproduce the 
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FIG. 4. Experimental (left) and simulated (right) FROG 
traces of the pulse propagating in the central waveguide mea- 
sured at the output of samples. Top: L = 25 mm. Center: 
L = 40 mm. Bottom: L = 60 mm. The wavelength scale is 
shown between the columns. Input power Po = 0.9 MW. 



observed features. The simulation reveals that the long- 
wavelength LB observed at z = 25 mm is red-shifted to 
Ac = 1950 nm (i.e. beyond the spectral acceptance of 
our FROG) while the bluest is actually the one observed 
at z = 40 mm. Further propagation to 2: = 60 mm leads 
to full decay of the LBs captured in the measurements 
at 2: = 25 and z = 40 mm, and only weaker, moderately 
red-shifted fragments surviving collisions in the central 
wavepacket are observed. 

The key to understand the LBs decay mechanism is to 
consider the influence of the wavelength-dependent dis- 
persion and diffraction parameters of our system. While 
propagating through the sample, the LB is red-shifted 
due to self-steepening and Raman scattering. In this 
course, the LB reshapes adiabatically to the increased 
dispersion and diffraction strength pertaining to the 
shifted wavelength. Because for longer wavelengths the 
coupling strength between waveguides increases exponen- 
tially, the energy threshold [7] required for the LB prop- 
agation grows steeply. Therefore the solitary wave will 
eventually decay when the dispersively modified energy 
threshold exceeds the energy of the excitation. The sce- 
nario is illustrated in Fig. 5. The curved surface defines 
the temporal width of the idealized LB solution of the dis- 
crete NLSE [12] without self-steepening and Raman scat- 
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FIG. 5. Green line: trajectory in parameter space of the 
brigthest LB appearing in the simulation of Fig. 3a. The 
blue-red surface shows the locus of the ideal stationary LB 
solution of the discrete NLSE. Energy /wavelength plane: re- 
gions of stability and existence of ideal LB; A-C indicate dif- 
ferent phases of the trajectory (see text for details). Gray 
lines: projection of the trajectory onto the orthogonal planes. 



point out that the identification of the decay mechanism 
opens new perspectives for the design of systems aimed 
at the excitation of very long-lived LBs. 

Definition of average diffraction and dispersion length. 
The average dispersion length over the path from zi to 
Z2 is defined as: 

II II 

where l3{oj) is the propagation constant of light in the 
single waveguide, ro{z) and 0Jc{^) are the 1/e pulse du- 
ration and the carrier angular frequency of the bullet at 
the propagation distance respectively. Analogously, we 
define the average diffraction length on the same path as: 

(i^diff) = / ^ dz (2) 

which relates to the waveguide coupling strength c{uj) 
pjj]. For the simulation of Fig. 3(a): (Ldisp) =2.6 mm, 
(^diff) = 12.7 mm. 



tering, determined uniquely for each carrier wavelength 
A and energy E of the central waveguide component of 
the LB. Additionally, the trajectory of the brightest pulse 
appearing in the dynamic simulation in Fig. 3(a) is over- 
laid to the plot (green line) along with its orthogonal 
projections (gray Hues). After strong pulse reshaping (re- 
gion A), the trajectory relaxes to the surface of the ideal 
LBs (region B), until the redshift pulls it behind the sur- 
face, where no LB can exist (region C). Because of the 
finite energy threshold for the existence of 3D discrete- 
continuous LBs \?\ this decay scenario is quite generic for 
experimental systems. Notice that the trajectory of the 
bullet in the E—X plane indicates that radiation losses are 
rather weak. To show that the decay mechanism is inde- 
pendent from radiation losses, we performed a simulation 
in an almost ideal Kerr medium with pure quadratic chro- 
matic dispersion, and longer coupling length (Ldiff(1550 
nm)=96 mm). These artificial conditions ensure that 
radiation from the LBs is strongly suppressed by large 
phase mismatch [17] and lower power excitation thresh- 
old, due to weaker coupling. The red shifting mechanism 
is retained by including self-steepening and Raman scat- 
tering. The results of the simulation reveal essentially 
the same dynamics as depicted in Fig. 5, while the prop- 
agation range of the LB normalized to Ldiff is expanded 
only by a factor of 1.5 (3 (i^diff) in total). 

In conclusion, we reported for the first time the ex- 
perimental observation of 3D nonlinear LBs. Our work 
revealed a new evolution scenario of adiabatic reshap- 
ing of LBs which sheds new light onpreviously puzzling 
results obtained in similar systems [l3|, and offers an in- 
terpretation in terms of evolving LBs for the dynamic 
of the so-called spatiotemporal 'quasi-solitons' [18j. We 
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